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ABSTRACT 

We report observations of gaseous methanol in an edge-on torus surrounding the young stellar 
object L1551 IRS5. The peaks in the torus are separated by ~ 10,000 AU from L1551 IRS5, and 
contain ~ 0.03 M s of cold CH3OH. We infer that the CH3OH abundance increases in the outer part 
of the torus, probably as a result of methanol evaporation from dust grain surfaces heated by the 
shock luminosity associated with the shocks associated with the jets of an externally located x-ray 
source. Any methanol released in such a cold environment will rapidly freeze again, spreading CH3OH 
throughout the circumbinary torus to nascent dust grains, planitesimals, and primitive bodies. These 
observations probe the initial chemical conditions of matter infalling onto the disk. 
Subject headings: ISM: general - ISM: individual (L1551, L1551 IRS5, HH 154) - ISM:jets and outflows 



1. INTRODUCTION 

In the earliest stages of their formation, low-mass 
young stellar objects are embedded in flattened gaseous 
envelopes. Surrounding these nascent protostars, plane- 
tary systems, primitive bodies and their attendant dust 
grains will condense from the infalling material. In this 
paper, we study the distribution of the organic molecule 
methanol around a protostellar core, as it is an important 
constituent of the young material in the disk. The L1551 
molecular core (distance 140 pc) contains a young binary 
system (seen as an infrared and radio source L1551 IRS- 
5) hidden inside a dense envelope providing ~ 150 mag- 
nitudes of visual extinction (Stocke et al. 1988, White et 
al. 2000, Fridlund et al. 2005). A molecular outflow em- 
anates from the center of the disk (Sncll et al. 1980, Kaifu 
et al. 1984, Fridlund & White 1989, Rainey et al. 1987, 
Parker et al. 1991), with atomic jet(s) associated with an 
x-ray emitting region being observed to the SW (Fridlund 
& Liseau 1998, Favata et al. 2002, 2003, Bally et al. 2003, 
Fridlund et al. 2005) . The protostellar disk is surrounded 
by a massive (radius ~ 20,000 AU) cool envelope that 
exhibits both rotational and infall motions (Takakuwa 
et al. 2004, Moriarty-Schieven et al. 2006). Fridlund et 
al. (2002) used HCO+, H 13 CO+ and 13 CO J= 1 - ob- 
servations to estimate that the mass of the disk is ~ 
2.5±1.5 M Q . 

2. OBSERVATIONS 

Observations of CH 3 OH (96 GHz 2 k -l k and 242 GHz 
5 fc -4 fe ), HCN (J= 1-0), and CO (J= 1 - 0) were ob- 
tained using the Onsala 20m and JCMT 15m telescopes 
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in Feb 2004 and Jan/ Apr 2005. The Onsala observations 
used an single sideband SIS receiver, with a 1600 chan- 
nel x 25 kHz autocorrelator backend. The JCMT obser- 
vations of a higher CH3OH transition (241 GHz 5fe-4/c) 
were made using the standard facility RxA2 receiver. All 
lines were corrected to a common main beam brightness, 
T mb scale. Additional data for CS and C 34 S ( J= 2 - 1 ) 
are taken from Fridlund et al. (2002). 

3. METHANOL 

Methanol plays a key part in the chemistry leading to 
the production of biogenic molecules, and is an impor- 
tant constituent of icy grain mantles and comets. Ther- 
mally excited gas phase CH3OH has previously been ob- 
served from warm inner envelopes surrounding just a few 
Class protostars, where it is evaporated from warm dust 
grains heated by the central protostar, or by an accretion 
shock at the edge of an infalling envelope (Goldsmith et 
al. 1999, Velusamy et al. 2002, Maret et al. 2005). 

Maps showing the distribution of CH3OH, HCN and 
CO toward L1551 observed from Onsala are shown in 
Figure 1, and velocity channel maps for the 96 GHz 
CH3OH and CS J= 2 - 1 lines in Figure 2. 

The methanol has three peaks located ~ 6500-12000 
AU from IRS5 along a SE-NW line. The high density 
tracer, HCN is centrally peaked and marginally resolved, 
located between the 2 inner methanol peaks. CO emis- 
sion, which traces low density material, streams orthog- 
onally away to the NE and SW. The methanol appears 
as an edge-on torus surrounding the HCN core, and is 
oriented orthogonally to the CO outflow. The velocity 
shift between the two CH3OH peaks is 0.07 km s _1 - 
which if indicative of rotation, would imply a rotational 
period ~ 6.5 x 10 5 yr. Spectra toward L1551 IRS-5 and 
at the two main methanol peaks are shown in Figure 3, 
along with the model fits described next. 

We have modeled the CH3OH line intensities using an 
accelerated lambda iteration radiative transfer technique 
for a spherical-core envelope (Phillips & Little 2000), 
with an assumed core radius of 5000 AU; envelope ra- 
dius of 12,500 AU; and turbulent velocity of 0.45 km s _1 
(from the H 13 CO+ linewidth - Fridlund et al. 2002). The 
A and E symmetry states of CH 3 OH were calculated sep- 
arately using collisional rates, energy levels and Einstein 
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Fig. 1.— CH 3 OH 2 -lo(A) line map (Onsala data at 96.7414 
GHz, green image with white contours) superimposed with con- 
tours of the HCN J= 1—0 (Onsala data shown as black contours), 
and (red) and (blue) shifted CO J= 1-0 outflows (shown as red and 
blue contour lines respectively) from the data shown in Fridlund 
et al. (2002). West is on the left, north at the top. The bottom left 
bar shows the size scale, for a distance of 140 pc. 



rates from Muller et al. (2001) and Pottage et al. (2004), 
and convolved to the Onsala and JCMT beamwidths. A 
number of models were run spanning a range of tem- 
perature, density and abundance gradients to simultane- 
ously match both the relative and absolute spectral line 
intensities, and the observed morphology of the torus. 
The observations unfortunately do not provide strong 
constraints on the CH3OH kinetic temperature, and so 
the temperature and density gradients used in the model 
were chosen to match those of White et al. 2000 derived 
from dust continuum measurements. 

Using this model the best fit suggests a central hydro- 
gen gas density (at the position of L1551 IRS-5) n(H 2 ) = 
7.5 x 10 5 cm -3 and a methanol abundance X = 8x 10~ 10 
relative to H2. Since the temperature and density are 
known to decline away from L1551 IRS5 (see Figure 4a, 
which is derived from White et al. 2000), the observed 
torus structure then requires that the methanol abun- 
dance increases to ~ 4 x 10 -9 at the inner edge of the 
methanol ring, and to 10~ 8 in the outer envelope, as 
shown in Figure 4b. 

The abundance increase is needed since H2 density is 
~ 20 times lower at the edge of the CH3OH peaks than 
towards the source centre at L1551 IRS-5. In this best 
fit model, the torus contains 0.03 of CH3OH. Mod- 
els run with CH3OH kinetic temperatures in the torus 
>20K failed to reproduce the observed disk structure 
whilst simultaneously matching the absolute and relative 
line temperatures. Since the gas and dust should be well 
thermalised at the estimated densities, it is reasonable 
to assume that the kinetic temperature of the CH3OH 
will be close to the dust temperature (13K from White 
et al. 2000). The models underestimate the observed 
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Fig. 2. — [upper] Channel maps in the 2§-\§{A) CH3OH line. 
The lowest contour level is 0.15 K km s — 1 and contour intervals 
are at 0.1 K km s _1 , [lower] Channel maps in the CS J= 2 — 1 line 
with the lowest contour at 0.375 K km s — 1 channel, and subsequent 
contour intervals increasing by 0.125 K km s _1 . The central (0,0) 
position of the maps and spectra is that of L1551 IRS5 at RA(2000) 
= 04 h 31 m 34 s . 13, Dec(2000) = 18° 08' 04". 95. 



linewidth toward IRS-5, which we believe is affected by 
additional dynamics and/or turbulence/outflow activity 
near L1551 IRS5 which is not modeled here. 

4. OTHER SPECIES IN THE CIRCUMBINARY TORUS 

The integrated HCN, CS, C 34 S, HCO+ and H 13 CO+ 
distributions peak close to L1551 IRS-5 (Fridlund et 
al. 2002), whereas the CH3OH map appears as an edge- 
on torus surrounding the HCN core. The CS and 
HCO + channel maps show a faint signature of this 
torus structure, but at low levels not contributing sig- 
nificantly to their integrated emission maps (Fridlund et 
al in preparation). 

Van der Tak et al. (2000) suggest that gas phase pro- 
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Fig. 3. — CH3OH data and model fits for the (left to right) 
2_i-l_i(E), 2 -l (A) and 2 -l (E) lines at 96.7394, 96.7414 and 
96.7446 GHz, and the 5_i-4_i(£), and 5 -4 (A), lines at 241.7672 
and 241.7914 GHz respectively toward the central IRS5 position, 
and at the two peaks in the 96 GHz CH3OH line. The frequency 
scales assume v; sr = 6.5 km s" 1 (Fridlund et al. 2002) 



duction of CH3OH at temperatures <100 K follows the 
radiative association of H2O and CHjj\ However, the rel- 
ative inefficiency of the process would result in low abun- 
dances ~ 10~ n relative to H2. Recent storage ring mea- 
surements (Gcppcrt ct al. 2006) of the dissociative recom- 
bination of CH 3 OH 2 + with electrons indicate extensive 
C-0 bond fragmentation, emphasising the difficulty for 
ion-molecule reactions to form CH3OH, which will prefer- 
entially be destroyed via an + CH3OH — > CH 3 OH^ 
+ c+ reaction. It is increasingly difficult to reconcile ob- 
served CH3OH abundances with ion- molecule chemical 
production. Solid phase reactions, including cosmic ray 
irradiation processes therefore likely are important for 
CH3OH formation on grain surfaces (Wada et al. 2006). 
Observations in dark clouds, where grain chemistry op- 
erates, indicate higher abundances of ~ 10 -9 (Takakuwa 
et al. 1998, Turner 1998), and toward hot cores, where 
CH3OH is evaporated off grain surfaces, abundances > 
10~ 7 (Blake et al. 1987) have been reported. Grain tem- 
peratures would need to be > 120 K for methanol to 
be efficiently released from icy or clathrate mixtures on 
dust grain surfaces (Blake et al. 1987, Sandford & Alla- 
mandola 1993). 

4.1. Heating of grains 

Our model suggests that gas phase methanol in the en- 
velope is sufficiently cold to rapidly freeze back onto dust 
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Fig. 4. — (a) Temperature and density from White et al. (2000) 
used in the CH3OH model, (b) abundance and column density 
estimated from the best fit CH3OH model 



grains and solid surfaces, which will occur on a timescale 
(t~3x 10 9 /n(iJ) - 4000 years, where n(H) is the 
hydrogen density cm -3 ). Shocks and direct heating, 
could re-evaporate the icy mantles, although the liber- 
ated CH3OH will rapidly cool again. The mid- plane dust 
temperature, Tdust due to L1551 IRS5 will be (Fridlund 
et al. 2002): 
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At the CH3OH peaks, T dust - 13 K (White et 
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al. 2000), which is inadequate for CH3OH to be evap- 
orated. Methanol has been observed toward the L1157 
molecular disk (Velusamy et al. 2002), and higher mass 
protostellar cores (van der Tak et al. 2000), where di- 
rect heating could directly evaporate the CH3OH - but 
it does not appear that this mechanism could drive off 
significant amounts of CH3OH at distance of the L1551 
CH3OH peaks. 

An alternative source of heating is suggested by the 
presence of an Hcrbig-Haro jet (HH154) emanating from 
L1551 IRS5, and terminating in a working surface - 
presently located ~ 15" away from the IRS5 binary 
(Fridlund et al.2005). This jet is one of the few X- 
ray sources associated with low mass stellar jets, and 
has 4 separate shock interactions - two detected in X- 
rays (Favata et al.2002, 2006, Bally et al. 2003), and 
two seen in [OIII] 5007A and Ha (Fridlund et al. 2005). 
Assuming a clear line of sight to the surface of the 
flared disk/envelope/torus surface, the temperatures of 
icy grains heated by the jet luminosity would be (Hollcn- 
bach & McKee, 1979): 



T - — 47 

J- grain ^ 1 



Juv 



K 



(2) 



for ice, where Ci is 2, the grain radius, a, n is assumed to 
be 0.15^m (Brown, 1990) and the mean intensity is Juv- 
The shock luminosity is then: 



F s = 5.8 x 1Q- A n a v 3 s7 



erg cm 



(3) 



The projected distance of the methanol peaks from the 
shocks is w 6000 AU. Assuming this flux is dominated by 
X-rays, EUV and FUV emission (Hollenbach & McKee 
1979; Hartigan et al. 1987) T gram will reach 100-150 K - 
which is adequate to liberate CH3OH. The grains would 
however rapidly cool again, with the evaporated gas be- 
coming thermally well coupled with the grains. The in- 
fall of material would naturally move grains initially at 
the edge of the torus from the zone exposed to x-rays 
and EUV/VUV photons to cooler shielded regions of the 
torus. 

The gas phase CH 3 OH abundance in L1551's outer 
molecular envelope will depend on the balance between 
several processes including: o) freeze-out onto dust 
grains; b) destruction by ion-molecule reactions, since 
enhanced ion densities will be induced by the x-ray irra- 
diation forming CH 3 OH on a similar timescale to grain 
surface reactions (Wada et al. 2006); and c) evaporation 
or release from dust grains following grain heating, or 
vapourisation by an accretion shock (Velusamy, Langcr 
& Goldsmith 2002). 

To test the importance of b) above, a chemical model 
was kindly run for us by Dr S. Viti (Lintott & Viti 2006 
in preparation) for a molecular cloud having solar metal- 
licity and a high (100 x standard) ionization rate, with 
the same temperature and density to those we infer. This 
suggests that a high methanol abundance could survive 
against ion-molecule reactions for at least 10 4 yrs i.e. 
similar to or greater than the freeze-out timescale. For 
c), we favor release from heated grains over the accre- 
tion shock scenario on the basis of the narrow CH3OH 
linewidths that are observed in the outer L1551 molecu- 
lar envelope (~ 0.6 km s _1 compared to several km s _1 in 
the L1157 inner core), and the distance of the CH 3 OH 



peaks from L1551 IRS-5. There is no compelling evidence 
to suggest the CH3OH abundance in the outer L1551 
torus is significantly influenced by an accretion shock: 
the narrow CH3OH line center lies within 0.1 km s _1 of 
the systemic velocity inferred from many other lines (~ 
6.5 km s" 1 - Fridlund et al. 2002). 

4.2. Discussion and Conclusions 

Star formation occurs in the dense cores of molecular 
clouds. Assuming that comet formation occurs as an 
adjunct to this process, information about the primor- 
dial chemical abundances should be locked into the the 
dust grains and comets, which are the main repositories 
of primitive material left over from the solar protostellar 
disc. The molecular inventory of protostellar envelope 
material has been discussed by White et al. (2003), with 
CH3OH commonly being found in comets (Mumma 
et al. 2003), along with other organic molecules. Our 
detection of methanol in a region at a similar distance 
from LI 551 IRS5 as that of the Oort Cloud from our 
own Sun, suggests that the CH 3 OH may be able to 
accrete to and spread throughout the outer L1551 
protostellar envelope, seeding dust grains and primitive 
bodies with chemical precursor material that could 
contribute to the synthesis of more complex molecules. 
This detection of methanol in the LI 551 circumbinary 
torus provides a first glimpse into the initial chemical 
conditions of matter infalling onto protostellar disks, 
and future determinations of the abundances of other 
molecules in the X-ray/EUV illuminated zone could 
therefore provide a interesting test of the "interstellar" 
versus "nebular" aspects of cometary and planetesimal 
chemistry. 

In summary: a) we have detected a massive (0.03 
M®), and most likely cold (< 20K - based on dust 
continuum measurements), toroidal shaped ring of 
methanol surrounding the L1551 IRS5 protobinary 
system, b) we propose a new mechanism where shock 
luminosity heating by the jets associated with an 
externally located x-ray source - these have a direct 
line of sight to the torus surface could raise the dust 
temperature sufficiently to liberate CH3OH, c) we 
suggest that the luminosity from the jet /x-ray source 
located just above the L1551 torus is able to heat the 
dust grains sufficiently to liberate a substantial amount 
of gaseous CH3OH in the outer disk of L1551, and d) we 
speculate that the methanol released will rapidly freeze 
back onto solid material in the outer envelope, spreading 
CH3OH to the surfaces of dust grains and primitive 
bodies - and potentially modifying the primitive surfaces 
of solid material in L1551's outer molecular envelope. 
The observations give indications of the initial chemical 
conditions of matter infalling onto the disks. 
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